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ABSTRACT 

HI in galaxies traces the fuel for future star formation and reveals the effects of feedback on neutral 
gas. Using a statistically uniform, H i-selected sample of 565 galaxies from the ALFALFA Ha survey, 
we explore HI properties as a function of star formation activity. ALFALFA Ha provides i?.-band and 
Ha imaging for a volume-limited subset of the 21-cm ALFALFA survey. We identify eight starbursts 
based on Ha equivalent width and six with enhanced star formation relative to the main sequence. 
Both starbursts and non-starbursts have similar Hi to stellar mass ratios (Mhi/M*), which suggests 
that feedback is not depleting the starbursts’ Hi. Consequently, the starbursts do have shorter Hi 
depletion times (tdepl implying more efficient Hi-to-H 2 conversion. While major mergers likely drive 
this enhanced efficiency in some starbursts, the lowest mass starbursts may experience periodic bursts, 
consistent with enhanced scatter in tdep at low M*. Two starbursts appear to be pre-coalescence 
mergers; their elevated Mhi/M* suggest that Hl-to-H 2 conversion is still ongoing at this stage. By 
comparing with the GASS sample, we find that tdep anti-correlates with stellar surface density for 
disks, while spheroids show no such trend. Among early-type galaxies, tdep does not correlate with 
bulge-to-disk ratio; instead, the gas distribution may determine the star formation efficiency. Finally, 
the weak connection between galaxies’ specific star formation rates and Mhi/M* contrasts with the 
well-known correlation between Mm/M* and color. We show that dust extinction can explain the 
Hi-color trend, which may arise from the relationship between M*, Mm, and metallicity. 


1. INTRODUCTION 

Gas inflows and outflows drive galaxy evolution by con¬ 
trolling the raw material from which stars form. The 
star formation history of the Universe may reflect the 
history of gas ac c retion onto dark matter halos (e.g., 
iKeres et all 120051 IProchaska &; Wolfe! 12009( 1. and ob¬ 
served galaxy scaling relations may likewise trace the 
history of gas flows. For instance, the observed rela¬ 
tions between galaxy stellar masses, star formation rates 
(SFRs), and metallicities may stem from variations in 
the e fficiency with wh i ch galaxies accrete and expel gas 
(e.g.. |Dalcantonll2007HMannucci et alJl20lfl iDave et alJ 

l2011bMal: ILillv et all 1201. llT~ 

The connection between gas content and star for¬ 
mation is a well-established result. Individual stars 
form from dense cores within molecul ar clouds (e.g., 
iMvers fc Bensonlfl983HMotte et al.lll998D . Onkiloparsec 
scales, a galaxy’s SFR surface density, Esfr, increases 
with the H 1 +H 2 gas surface density, Shtlh , , as pa¬ 
ra mete rized by the Kennicutt-Schmidt Law dKennicuttl 
Il998bi) . More recent work shows that £h 2 , rather 
than Shi, drives the Kennicutt-Schmidt la w, even in 
Hi-dominated regimes (jSchruba et al.l[20lTh . Above a 
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threshold density of ~10 M 0 /pc 1 2 , Hi “saturates”; this 
column density is sufficient to shield molecular gas from 
photodissociatio n. Most gas above this threshold i s 
molecular (e.g., iWong fc B1 it zl [20021 : iBigiel et al .1120081) . 
resulting in no trend between SFR and HI in this den¬ 
sity regime. Below this threshold, however, Esfr and 
Xhi correlate, albeit with a l arge scatter , d ue to the re¬ 
lation between Shi and Sh 2 dSchruba et al.l[20TTh . 

Compared to the link between H 2 content and SFR, 
the relationship between galaxies’ HI content and star 
formation is less straightforward. Galaxies’ Hi typically 
extends to much l arger radii than the stellar distribu¬ 
tion (e.g., iBroeils fc Rheel Il997f ) and may constitute 
a gas reservoir for fueling future star formation. Ac¬ 
cretion of gas from the intergalactic medium may re¬ 
plenish this reservoir, and gas flows may bring HI in¬ 
ward, lea ding to star form ation in the inner regions of 
galaxies. IProchaska fc Wolfe! ( 2009 ) suggest that galax¬ 
ies’ Hi disks exist at a constant, unstable density, with 
any subsequent accretion leading to the creation of stars 
and resulting in an SFR that traces the accretion rate. 
In massive galaxies, the total HI gas fraction correlates 
with signs of recent accretion, such as an outer metallicity 
drop, and this accretion appe ars t o p ower s tar formation 
throughout the galaxy disk dMoran et al.l 1201 2). The 
HI in low-mass galaxies, on the other hand, may not sig¬ 
nify recent accretion. Low-mass galaxies tend to be more 
gas-rich than high-mass galaxies, and HI constitutes the 
dominant component of their interstellar medium (ISM). 
The long gas consumption times in such systems may in¬ 
dicate that_starformation proceeds inefficiently. Alterna¬ 
tively, [Kannappanejnil] (I2Q13T) argue that star formation 
cannot keep pace with the rate of gas accretion, causing 
HI to accumulate. The relationship between HI and star 
formation may also change in dwarf galaxies due to their 
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lower metallicities. At low metallicities and consequently, 
low dust content, the format ion of a given Yh-, re quires 
a higher Hi column density (|Krumholz et al.112009T) . As 
a result, the relationship between Hi density and H 2 for¬ 
mation differs for low-metallicity galaxies, with galaxies 
such as the Small Magellanic Cloud exhibiting a higher 
threshold density for Hi saturation dBolatto et al.ll201lh . 

The advent of large Hi surveys capable of resolv¬ 
ing individual galaxies has clarified the connection 
between H 1 mass and galaxy properties. In par¬ 
ticular, the Arecibo Legacy Fast ALFA (ALFALFA) 
survey is a blind 21-cm survey, which covers 7000 
deg 2 and has detected ~ 3 0,000 galaxies out t o z = 0.06 
(|Giovanelli et all[2005bl ial: lHavnes et all 1201 111 . In con¬ 
junction with_^ALFALFA 1 _the GALEX Arecibo SPSS 
(GAS S; ICatinella et al.112010H and Ha3 (IGavazzi et al.1 
I2012H surveys have investigated HI and star formation in 
different galaxy regimes. The GASS survey examines the 
Hi content of massive galaxies (M* > 10 10 M 0 ) using 
additional Arecibo observations for galaxies undetected 
in ALFALFA. Ha3 studies the effect of environment on 
Hi content and star formation in ALFALFA galaxies 
within the Local Supercluster. The ALFALFA survey 
and related surveys have established scaling relations be¬ 
tween H 1 gas fraction and galaxy stellar mass, stellar sur¬ 
face density, color, SFR, and specific SFR (sSFR; e.g., 
Catinella et al.1120101 iHuang et al.1 l2012t IGavazzi et al'l 
2013 1. The positive correlations found between Hi gas 
fraction and blue color or sSFR imply a link between 
galaxies’ HI content and their current global star forma¬ 
tion. 

Starburst galaxies may depart from the typical rela¬ 
tions between HI content and star formation, however. 
The Kennicutt-Schmidt Law may differ for starburst 
galaxies, with starburst galaxies forming stars m ore effi¬ 
ciently from a given mole c ular gas mass (e.g., iKennirnittl 
Il998bt iDaddi et all 120101 ; iGenzel et al.ll2010f) . One pos¬ 
sible explanatio n for this increase d efficiency is merger 
activity (e.g., lYoung et all 119861 : iSanders et all 119861 : 
iGombes et all 11994( 1 . and many starbursts appear to be 
interacting systems. If starbursts gain H 1 gas via major 
mergers instead of accretion, they may differ from non- 
starbursts in both their H 1 gas consumption times and 
Hi gas fractions. In addition, due to their young stel¬ 
lar populations, mechanical and radiative feedback will 
have a stronger effect on the ISM of starbursts. This 
feedback may decrease the H 1 gas fractions of starbursts 
Iw driving outflows or ionizing the neutral gas. lOev et ahl 
(I2007T) suggest the latter scenario as an explanation for 
the lower H 1 gas fractions in starburst galaxies in the 
Survey for Ionization in Neutral Gas Galaxies (SINGG). 
The SINGG result contrasts with the ALFALFA and 
GASS trends of higher HI gas fractions in more highly 
star forming galaxies and demonstrates that the H 1 con¬ 
tent of starbursts requires further study. 

Starburst galaxies present an opportunity to study the 
relations between neutral gas content, star formation, 
and feedback in extreme conditions. The Hi gas frac¬ 
tions and kinematics of starbursts may reveal the mech¬ 
anisms for triggering extreme star formation episodes 
and the impact of feedback on global gas content. Pre¬ 
vious studies of H 1 in starbursts have focused on in¬ 
dividual galaxies or optically selected samples (e.g. , 
lYun et al.lll993t iHuchtmeier et al.l[2007t lOev et al.ll2007l: 


iLonez-Sancliez et all l2012ll . To systematically compare 
the H 1 properties and star formation efficiencies of gas- 
rich starbursts and gas-rich non-starbursts, we use the 
“Fall-sky” portion of the ALFALFA Ha survey. AL¬ 
FALFA Ha is a volume-limited subset of the ALFALFA 
survey consisting of 1555 galaxies with follo w-up Ha and 
R-band imaging (lYan Sistine et al.1 l2015h . With the 
565 galaxies in the completed “Fall sample” of the AL¬ 
FALFA Ha dataset, we investigate the regulation of the 
H 1 gas supply throughout the star-formation process. 

2. DATA AND METHODS 
2.1. The ALFALFA Ha Survey 

The recently-completed ALFALFA survey is a blind 
21-cm survey with Arecibo that covers 7000 deg 2 of 
sky. The Hi-selected ALFALFA Ha survey consists 
of all ALFALFA-detected galaxies within two desig¬ 
nated areas, a Fall-sky regi on and a Spring-sky re¬ 
gion (lYan Sistine et al.fl2015H . The Fall sample only 
includes galaxies with reliable Hi detections (i.e., AL¬ 
FALFA codes 1 and 2) and with recession velocities 
v = 1460 — 7600 kms” 1 . These velocities correspond 
to distances of ~20-100 Mpc, and the sample is volume- 
limited for Mhi > 10 9 3 Mq. Hi masses, 21-cm velocity 
widths, a nd distanc es come from the ALFALFA catalog 
(jHavnes et al.ll2f)ITI) . The selection of the most p r obabl e 
optical counterparts is described in lHavnes et al.1 (|2011h : 
ambiguous optical identifications or blended H 1 signals 
from multiple galaxies occur approximately 10% of the 
time. 

In this work, we consider the complete, AL¬ 
FALFA Ha Fall sample, whi ch contains 565 galaxies ; 
the full sample is described in lYan Sistine et al.1 (I2015H . 
To compare the global star formation of these galax¬ 
ies with their global Hi content, we obtained .R-band 
and Ha imaging for the Fall-sky ALFALFA Ha galax¬ 
ies with the WIYN 0.9m telescope^ and the Kitt Peak 
National Observatory (KPNO) 2.1m telescope between 
Sept. 2006 and Oct. 2012. The parameters of the AL¬ 
FALFA Ha survey, observations, an d data reduction are 
described in lVan Sistine et all (|2015fl . At the 20-100 Mpc 
distances of our sample, the observed Ha emission typ¬ 
ically extends over several tens of arcsec in each galaxy, 
and the 3" spectroscopic fiber of the Sloan Digital Sky 
Survey (SDSS) does not accurately capture the total star 
formation. 

Ha and R-band fluxes were measured for each galaxy 
individually using aperture photometry. The Ha images 
were continuum-subtracted prior to flux measurement, 
and R-band fluxes were corrected for Ha contamination 
within the bandpass. Of the 565 Fall sample galaxies, 
542 galaxies are detected in Ha. To facilitate compari¬ 
son with the literature, the quoted Ha equivalent widths 
are an observed quantity and are not corrected for ex¬ 
tinction or [N 11 ] emission. For derived parameters, such 
as SFRs, the Ha fluxe s were first corrected f or Ga lac- 
tic absorption using t helSchlaflw fc Finkbeiner; (20J I: re¬ 
calibration of the lSchlegel et all (|1998H extinction maps. 

6 The WIYN 0.9m telescope is operated by WIYN Inc. on behalf 
of a Consortium of partner Universities and Organizations (see 
www.noao.edu/0.9m for a list of the current partners). WIYN 
is a joint partnership of the University of Wisconsin at Madison, 
Indiana University, the University of Missouri, and the National 
Optical Astronomical Observatory. 
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Ha fluxes were then corrected for [N n] contamination 
and internal absorption using the observed R-band abso¬ 
lute magnitudes (Mr) and scaling relations derived from 
a sample of 803 star-forming galaxies from_the_ K PNO 
Internation al Spectroscopic Survey (KISS; ISalzer et al.l 
200011200511. The s e corr ections are described further in 
Van Sistine et al.l (l2015h . We compare these extinction 
corrections to corrections derived from Wide-fie ld In¬ 
frared Survey Expl orer ( WIS1 IWright et al.ll20ldT) pho¬ 
tometry in § HI After correcting the Ha emission 
for Galactic and internal extinction and [N ii] emission, 
we convert the Ha luminosit ies to SFRs. The AL ¬ 
FALFA Ha SFRs present ed in lYan Sistine et all (|2015l) 
use the iKennicuttl (11998al) calibration, which assumes a 
ISalpeteiT l 19551) initial mass function (IMF). Here, we 
scale the [Kennicutt (1998a) calibration to a iChabrieil 
(|2003D IMF and calculate the SFR as: 

SFR = 4.6 x 10~ 42 L(Ha), (1) 

where L(Ha) is the Ha luminosity in ergss -1 and the 
SFR has units of M 0 yr _1 . The iChabrierl (120031) IMF 
results in SFRs th at are a factor of 1.7 lower than SFRs 
calculated with a iSalpeterl (119551 ) IMF. 

2 . 2 . Stellar Mass Estimation 

To estimate stellar masses for our sample, we 
need to account for galaxy-to-galax y variations in 
the st ellar mass-to-light ratio ( M/L ). iBell &; de Jond 
(I2001D demonstrate that galaxy M/L ratios should vary 
systematically with star formation history and enrich¬ 
ment history, leading to a dependence of M/L on galaxy 
color. Accounting for the correlation between M/L ratio 
and color is particularly important to accurately esti¬ 
mate the masses of starburst galaxies, whose young stel¬ 
lar populations amplify the luminosities of all UV-NIR 
bands. 

To correct for this effect, we o btain g ala xy colo rs from 
the SDSJ0 Ninth Data Release (lAhn et al.ll2012fl . SDSS 
data are available for 513 of the 565 galaxies in the 
Fall sample. However, SDSS photometry can be prob¬ 
lematic, particularly for low surface brightness or ir¬ 
regular galaxies. The SDSS de-blending pipeline sep¬ 
arates overlapping objects and sometimes incorrectly 
shreds one “parent” galaxy into multiple “ children” 
(e.g., lAbazaiian et al.lT2004l IWest et aLll 2010 D . To de¬ 
termine whether shredding is a concern for our sample, 
we examine the SDSS data by eye for 10% of the AL¬ 
FALFA Ha galaxies. For 90% of these galaxies, the g 
and r photometry of the brightest de-blended child and 
the parent object agree to within 0.2 mag. We there¬ 
fore select the brightest de-blended child for each AL¬ 
FALFA Ha galaxy from the SDSS catalog. As an addi¬ 
tional check, for the full ALFALFA Ha Fall sample, we 
compare the SDSS r-band magnitudes with our R-band 

7 This publication makes use of data products from the Wide- 
field Infrared Survey Explorer, which is a joint project of the Uni¬ 
versity of California, Los Angeles, and the Jet Propulsion Labo¬ 
ratory/California Institute of Technology, funded by the National 
Aeronautics and Space Administration. 

8 Funding for SDSS-III has been provided by the Alfred P. Sloan 
Foundation, the Participating Institutions, the National Science 
Foundation, and the U.S. Department of Energy Office of Science. 
SDSS-III is managed by the Astrophysical Research Consortium 
for the Participating Institutions of the SDSS-III Collaboration. 


photometry. The photometric data from SDSS and the 
ALFALFA Ha R-band data show a tight, linear relation¬ 
ship (Figure [T]), although noticeable outliers exist. An 
examination of the outliers shows that they are caused 
by SDSS de-blending errors, overlapping or nearby bright 
stars, and incorrectly separated galaxy pairs. We treat 
the SDSS photometry as unreliable if it differs from a 
least-squares fit to the R-band data by more than 1 mag 
(Figure [I]). We also eliminate the SDSS photometry for 
one additional galaxy with a discrepant 5 -band magni¬ 
tude that results in an unrealistically red g — r color. We 
do not include these galaxies in any analyses that rely 
on SDSS-derived parameters, such as stellar masses or 
galaxy radii. These cuts lead to a sample of 489 AL¬ 
FALFA Ha galaxies with SDSS photometry. Including 
the galaxies with uncertain SDSS photometry does not 
affect any of our conclusions in the following sections. 

We calculate stellar masses using the SDSS r-band 
luminosities and the relationship between M/L ratio 
and g — r color from Table 7 of IBell et al.l (120031) . We 
first subtract the observed Ha fluxes from the r-band 
fluxes to ensure that nebular emission does not affect 
the galaxy colors or magnit udes. In addition, we sub¬ 
tract 0.093 dex from the IBell et ahi (120031) M/L ra¬ 
tios to convert from a “di et” iSalpeted (119551) IMF (see 
IBell fc de Jond [20011) to a IChabrieil (l2003D~Initia l mass 


function ( Gallazzi et al.1 [20081 : IZibetti et all [2009D . Sys¬ 
tematic uncertainties from the assumed dust extinction 
and s tar formation his tories dominate the un certainties 
in thelBell et, al.1 (I2003D MIL ratios. Following lBell et al.1 
(I2003D . we adopt a total systematic uncertainty of 0.1 dex 
for the M/L estimates. We also compare the resulting 
st ellar mass est i mates with th e alternative pre scription 
of IZibetti et al.1 (|2009). While IBell et al.l (|2003f ) assume 
a s mooth star formation history to model galaxy colors, 
IZibetti et al.l (|2009D consider th e effect of bu r sts. A s a 
result, masses obtained following I Zibetti et al.l (|2009D are 
generally lower than the IBell et alT ( 20031 ) mass estimates 
(Figure [2]). In particular, at the low-mass end, the stellar 
masses differ by about a factor of three. However, the 
relative stellar masses of the gal axies in th e sa mple are 
only weakly affected. We use the lBell et al.1 (2 003i ) stellar 
mass estimates for th e rest of our analysis, but we note 
that using the IZibetti et al.l (120091) mass estimates does 
not change our conclusions. 


2.3. Selection of Starbursts 

The existing literature contains a variety of defini¬ 
tions as to what cons t itutes a starburst g alaxy (e.g., 
iKnapen fe Jamesl 120091 : iBergvall et al.l 120151) . Common 
starburst definitions are based on high sSFRs or Ha 
equivalent widths (EWs), high sSFRs relative to simi¬ 
lar mass galaxies, high SFRs per unit area, or short gas 
consumption times. In addition, some definitions only 
consider global star formation enhancements, while oth¬ 
ers consider smaller-scale enhancements as well. A nu¬ 
clear starburst, for instance, will have a high SFR in the 
nuclear region but may not necessarily qualify as a star- 
burst based on the galaxy’s total SFR. 

Low-redshift studies often define starbursts as galaxies 
with high SFRs relative to their past average rates (e.g., 
Ostlin et al.1 120011: Brin chmann et al.1 12004 iLee et akl 
20091 : IBergvall et aD 20151) . Here, we adopt a similar def- 
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Comparison of SDSS and ALFALFA Photometry 



Fig. 1.— A comparison of the SDSS r-band and the ALFALFA Ha .R-band photometry for galaxies in the Fall sample. The line represents 
a least-squares fit to all data points. Statistical errors are typically smaller than the symbol size. Asterisks indicate galaxies with r-band 
magnitudes that differ from the fitted line by more than 1 mag. 


Comparison of Stellar Mass Estimation Methods 



Fig. 2.— A comparison of the stellar masses estimated following IBell et all (!200. > i and IZibetti et aTli 11 20091 ). The solid line indicates 
a one-to-one relation. All masses were estimated using r-band lumi nosities an d g — r col ors. T he erro r s show n include the statistical 
photometric errors and the 0.1 dex scatter in the M/L -color relation (IBell et al.ll200 3L The IZibetti et all 12009 *) stellar masses are lower, 
particularly at the low-mass end, due to their adoption of bursty star formation histories. 


inition and define starbursts as galaxies with Her equiv¬ 
alent widths (EWs) greater than 80 A. This EW corre¬ 
sponds to a birth rate parameter of ~2, i.e., the star- 
bursts have instantaneous SFRs greater t han or equal to 
twice their past average SFR dLee et al.1 12009H . Since 
dust in starburst galaxies may attenuate the ionizing 
continuum radiation more than the optical stellar con¬ 
tinuum, dust effect s may lower the observed Ha EW s 
in starbursts (e.g., iCalzettilIl997fc ICharlot fe Kalll 120001) . 
The Ha EWs of the starbursts may therefore underesti¬ 
mate their true ratios of current to past star formation. 
Eight galaxies (1.4% of the sample) have EWs above the 
80 A cut (Figure [3]), which is equivalent to an sSFR cut 
of approximately 6 x 10 -10 yr -1 . We list the EWs, SFRs, 
and HI and stellar masses of the starbursts in Table [L] 


We also consider an alternative definition of starbursts 
as galaxies that fall above the main sequence. Figure [4] 
shows a least-squares fit line to the SFRs and stellar 
masses of the sample. An additional six galaxies have 
SFRs more than 2<r above the best-fit line, although none 
are strong outliers above the main sequence. Table [T] 
summarizes the properties of these galaxies. Hereafter, 
we refer to these six galaxies as the “High sSFR Star- 
bursts” to distinguish them from the eight “High EW 
Starbursts” identified based on Ha EW. 

Since we are incomplete below Mhi= 1O 9 ’ 3 M 0 , we 
may not detect all the starbursts in the Fall-sky volume. 
For a rough estimate of our completene ss, we examin e 
the la rger, optically-selected sample of Bothwcl l et al.1 
( 2009 1, which contains 1110 galaxies at distances less 
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TABLE 1 

Properties of the ALFALFA Ha Starbursts 


ID 

RA”'" 

Dec a ’° 

Distance** 

(Mpc) 

Ha EW 

(A) 

SFR 

( Mgyr -1 ) 

Log (Mhi) & 

(m q ) 

Log (M„)° 
(M e ) 

^dep 

(Gyr) 

A d 

W 20 /W° 0 

High EW 

Starbursts 

AGC 112546 

01:31:20.6 

+28:48:29 

66.4 

286.6+4.8 

0.49+0.23 

8.7+0.1 

7.2+0.1 

1.1+0.6 

0.24 

<1.32 

AGC 330517 

23:32:04.7 

+28:57:21 

78.4 

276.4+4.8 

6.2+3.0 

9.6+0.1 

8.0+0.1 

0.7+0.4 

0.68 

1.43+0.43 

AGC 122866 

02:11:31.4 

+24:12:51 

37.3 

124.5+4.3 

0.02+0.01 

8.3+0.1 

7.6+0.1 

8.4+4.4 


<1.54 

AGC 120193 

02:22:55.0 

+25:18:53 

62.7 

101.9+2.0 

0.63+0.30 

9.6+0.1 


7.1+3.7 

0.18 

1.24+0.16 

AGC 333529 

23:22:39.8 

+28:57:18 

85.2 

92.0+10.0 

0.08+0.04 

9.0+0.1 

7.8+0.1 

12.8+6.8 


1.69+0.28 

AGC 330500 

23:30:10.0 

+25:32:01 

82.2 

81.6+1.3 

3.9+2.0 

9.9+0.1 

9.3+0.1 

2.1 + 1.1 

0.43 

1.25+0.15 

AGC 122187 

02:09:55.7 

+27:32:25 

67.6 

81.3+6.1 

0.05+0.02 

9.0+0.1 

7.9+0.1 

19.3+10.1 


1.39+0.15 

AGC 331191 

23:28:48.9 

+24:52:10 

71.2 

80.4+2.3 

0.65+0.31 

9.6+0.1 

9.0+0.1 

5.9+3.1 

0.18 

1.62+0.33 


High sSFR 
Starbursts 


AGC 122420 

02:38:50.5 

+27:21:59 

19.5 

77.2+2.1 

0.03+0.02 

8.0+0.1 

7.3+0.1 

3.2+1.8 


< 1.33 

AGC 320466 

22:57:20.7 

+27:58:52 

43.3 

72.7+8.1 

0.06+0.03 

9.1+0.1 

7.7+0.1 

23.2+12.2 


<1.19 

AGC 102643 

00:21:36.7 

+25:28:58 

95.8 

66.2+1.3 

0.66+0.32 

9.3+0.1 

9.1+0.1 

2.8+1.5 

0.29 

1.13+0.15 

AGC 330186 

23:17:17.7 

+28:36:03 

96.3 

23.3+0.4 

1.97+0.98 

9.6+0.1 

9.7+0.1 

2.1 + 1.1 


<1.12 

UGC 470 

00:44:14.4 

+26:50:35 

72.7 

66.9+1.3 

2.58+1.24 

10.2+0.1 

9.8+0.1 

5.9+3.1 

0.29 

1.09+0.01 

UGC 12821 

23:52:23.6 

+28:46:15 

91.7 

40.7+0.7 

10.09+4.98 

9.9+0.1 

10.5+0.1 

0.7+0.4 

0.34 

1.40+0.10 


“Coordinates of the detected Hi sou rce. _ 

b Values from the ALFALFA catalog (jHavnes et al.ll201lTh 

c The errors listed include the statistical photometric errors and the 0.1 dex scatter in the M/L -color relation. 
“i?-band asymmetry. 


than ~43 Mpc. While the iBothwell et all (120091) sample 
selection is less uniform than the ALFALFA Ha sample, 
it does contain a higher fraction of low-mass galaxies. As¬ 
suming the lowest ra tios o f Mht/SF R obse rved for late- 
type galaxies in the IBothwell et al.l ( 2009 ) sample and 
converting to a ICliabrierl ( 20031 ) IMF, we should detect 
all starbursts with M* = 10 s M 0 and sSFR> 10~ 7 yr -1 
and all starbursts with M* = 10 9 M 0 and sSFR> 10 -8 
yr -1 . However, we may miss any low-mass starburst 
galaxies whose gas is predominantly molecular. Al¬ 
though the few low-mass starbursts with CO measure¬ 
ments generally a ppear to have larger Hi masses than 
H? masses (e.g., iKobulnickv et al. l 119951: IStil fo Israeli 

2002: iBravo- Alfaro et al.ll2004lLsrae]il2005tlNidever et al.l 

20131 ). we caution that we may not detect the most ex¬ 

treme, low-mass starbursts with the highest H 2 /H 1 ra¬ 
tios. 

2.4. WISE Data 

Dust extinction may substantially affect the observed 
star formation properties of the ALFALFA Ha galax¬ 
ies. Using infrared measurements from WISE , we verify 
the accuracy of the ALFALFA Ha extinction corrections. 
Later, in § 13.1.21 we use the WISE data to examine the 
relationship between dust and Hi content. 

We obtain 3.4, 4.6, 12, and 22 /im fluxes for the 
sam ple from the WISE All-Sky Release So urce Cata¬ 
log (jWright, et al.l 120101 : iCutri fc et al.1 l2012fl . Elliptical 
aperture photometry from WISE is available for sources 
identified in the Two Micron All Sky Survey (2MASS) 
Extended Source catalog dJarrett et al.1 12000|) . We use 
this elliptical aperture photometry, where it exists, for 
objects flagged as resolved sources in the WISE cata¬ 
log and the profile fit magnitudes for all other sources. 
We then match each ALFALFA Ha galaxy to its closest 
WISE co unterpart within 6", the resolution at 3.4 /tm. 
Following iJarrett et al.l (|2013l ) , we apply the necessary 
color and magnitude corrections to the WISE photome¬ 
try. Of the 565 ALFALFA Ha Fall galaxies, 263 galaxies 


have detections in all four WISE bands (3.4, 4.6, 12, and 
22 /rm) and an additional 95 galaxies have detections in 
three WISE bands (3.4, 4.6, and 12 /im). 

The extinction corrections adopted for the AL¬ 
FALFA Ha galaxies in § 12.11 are based on the galaxies’ 
l?-band luminosities. Since this is a statistical correction, 
it may not be correct on an individual galaxy basis. To 
test the adopted extinction correction, we use the WISE 
and Ha observations to esti mate ex tinctio n-co rrected 
SFRs for individual galaxies. IWen et ah ! (I2014D cali¬ 
brate the WISE 12 /im and 22 /rm bands as extinction 
indicators using the Ha/H/3 ratios from SP SS spectra 
o f z < 0.25 star-forming galaxies. Following IWen et al.1 
mb, we calculate 

SFR = 0.87 x 10- 41 ' 27 L H a,ob s + avL„, (2) 

where LHa.obs is the Ha luminosity without an inter¬ 
nal extinction correction, vL v is the appropriate WISE 
luminosity. The coefficient a depends on the WISE 
band and extinction law adopted and ranges from ~0.02- 
0.04. T he factor o f (187 in the equation converts from 
the IKennicutt fe Evansl (120121) SFR calibrat ion, which 
uses a IKroupa &; Weidnerl (1200311 IMF to the iKennicuttl 
(|1998aJl calibration and the iChabrieil (1200311 IMF used 
in §0 We calculate extinction-corrected SFRs using 
t he WI SE 12 urn and 22 /rm luminosities and assuming 
a lCalzetti et al.l (120001) extinction law. 

We compare these SFRs with the ALFALFA Ha SFRs 
in Figure [5] Overall, the WISE -derived SFRs and AL¬ 
FALFA Ha SFRs show reasonable agreement. A few 
galaxies show higher SFRs using the WISE extinction 
correction than the original ALFALFA Ha extinction 
correction. These galaxies show prominent dust lanes 
or appear reddened in SDSS images. For these partic¬ 
ular objects, the original extinction correction used to 
derive the ALFALFA Ha SFRs is likely insufficient. The 
WXS'E-derived SFRs appear systematically lower than 
the ALFALFA Ha SFRs by ~0.1 dex for many galax¬ 
ies, especially point-source WISE detections at moderate 
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Fig. 3.— The distribution of the Ha EWs in the ALFALFA Ha Fall sample. The solid region of the histogram indicates the identified 
starburst galaxies. The hatched bin at negative EWs represents galaxies without Ha detections. The inset shows a zoomed-in view of the 
high-EW section of the histogram. 
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Fig. 4.— SFR and M* for the ALFALFA Ha Fall sample. Color shows Mhi/M*, and diamonds indicate starburst galaxies, as identified 
by Ha EW. The black line shows a least-squares fit to the data. Triangles indicate additional starburst candidates, whose SFRs are at 
least 2cr above the best-fit line. The black cross at the lower right shows representative error bars. At a given M*, galaxies with higher 
SFRs tend to have slightly higher HI gas fractions. 
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SFRs. The WISE profile-fit fluxes for these sources may 
be underestimated, as discussed in lCutri et al.1 (I2012H . In 
general, however, the WISE-ba,sed extinction corrections 
agree well with the adopted ALFALFA Ha extinction 
corrections. 


3. RESULTS 

3.1. Hi Gas Content 

3.1.1. The Hi Gas Supply of Starbursts 

Star formation is intimately linked with the cold gas 
content of galaxies. A large Hi supply may be necessary 
to fuel high SFRs, but the resulting feedback may expel 
or ionize much of the Hi gas. Given the extreme levels 
of star formation and feedback in starburst galaxies, we 
consider whether their HI content differs from the other 
galaxies in the ALFALFA Ha sample. High Hi gas frac¬ 
tions in starbursts may suggest that a large gas reservoir 
is a key precondition for triggering a starburst, while Hi- 
deficiencies might indicate that radiative feedback plays 
the dominant role in shaping starbursts’ ISM. 

In non-starbursts, larger Hi reservoirs do appear 
to lea d to e nhanc ed star formation. As shown by 
IHuang et alJ (|2012l i for ALFALFA galaxies, SFR corre¬ 
lates with Mhi over almost 4 dex in Hi mass, indicating 
a link between atomic gas and star formation. Since part 
of this trend may result from the tight relation between 
SFR and galaxy stellar mass, we examine whether highly 
star-forming galaxies have more Hi than galaxies of a 
similar mass. In Figure U we show the ratio of Mhi to 
stellar mass (M*) as a function of SFR an d M* . For con¬ 
sistency with previous F t I studies (e.g., iCatinella et al.l 
l20lii IHuang et all [Ml . we refer to Mhi/M* as the 
Hi gas fraction. Previous studies have found that at a 
given stellar mass, gala xies with higher SFRs tend to 
be m ore Hi-rich (e.g., iWang et al.ll201ll : IHuang et all 
l20H . and the ALFALFA Ha sample likewise exhibits 
this trend. The larger , 40% complete ALFALFA sam¬ 
ple (IHuang et a1.ll20 12ll shows that this relation between 
SFR and HI content also holds at lower HI gas fractions 
than probed by the ALFALFA Ha galaxies. However, 
while some of the starbursts’ SFRs are an order of mag¬ 
nitude higher than other galaxies of the same mass, their 
HI gas fractions do not show a comparable increase. In 
fact, the starburst Hi gas fractions are similar to those 
of star-forming galaxies with much lower SFRs. In the 
starburst regime, Hi content and star formation do not 
appear to be closely coupled. 

The starbursts are Hi-rich relative to the AL¬ 
FALFA Ha sample as a whole, but this Hi-richness re¬ 
sults from their lower-than-average stellar masses. Low- 
mass galaxies in general tend to be mo re Hi-rich (e.g., 
iGavazzi et all Il99fit IHuang et all l2012Tl . and in most 
cases, the Hi gas fractions of the starbursts are similar 
to other galaxies of the same mass. The relatively high 
Mhi detection limit of the ALFALFA Ha survey leads 
to a tight apparent relation between Mhi/M* and M*; 
Figure [6] displays the least-squares fit to this relation for 
the non-starbursts in the sample. Most of the starbursts 
have Mhi/M* values that are within ltr of the best-fit 
line. Only two of the high EW starbursts and one high 
sSFR starburst have Mhi/M* ratios that appear high for 
their stellar mass, deviating by > 2tr from the main trend 
(FigurejG]). Two starbursts even have lower-than-average 


HI gas fractions. This fact suggests that the high SFRs 
of the starbursts are not caused by an excess of Hi, but 
rather by an enhanced efficiency of converting HI into 
molecular gas and stars. 

Although the starbursts are not excessively H i-rich for 
their stellar masses, they are also not H i-deficient. Thi s 
conclusion contrasts with the results of lOev et alJ (120071! 
for sta rbursts in the SINGG sample. However. lOev et all 
lj2007l) use R-band magnitude as a proxy for stellar mass 
and identify starbursts based on Ha surface brightness, 
rather than EW. To compare with the SINGG results, 
we calculate the SFR surface density for a flat disk mor¬ 
phology as 


Xsfr = 


SFRsq 

2nRl 0 ’ 


(3) 


where R 50 is the R-band half-light radius and SFR 50 is 
the SFR within that radius. We show the relation be¬ 
tween Xsfr and the ratio of Mhi to R-band luminos¬ 
ity, Mhi/Lr, in Figure [7jr. Our sample spans the same 
range of M m /L R as the SINGG sample, and we like¬ 
wise see that galaxies with higher Xsfr do appear to 
have lower ratios of Mhi/Lr. On the other hand, as 
discussed in § m R-band luminosities systematically 
overestimate the stellar masses of starbursts, leading to 
underestimates of their HI gas fractions. In Figure [Tjr 
we replace L R with M*, calculated as described in § 12.21 
The trend between Hi content and Xsfr now appears 
weaker, and even at the highest Xsfr end, galaxies span 
a range of Hi gas fractions. This result highlights the 
importance of multicolor photometry for accurate stel¬ 
lar masses and gas fractions. The weak residual neg¬ 
ative correlation between Mhi/M* and Xsfr in Fig¬ 
ure is due to the scalings of Hi gas fraction, SFR, 
and half-light radius with M* in the ALFALFA Ha mass 
regime ('e.g.. ICatinella et alJ 120101: IHuang et alJl2012t 
iBrinchmann et alJ 120041 : iShen et al l l2003fb The large 
range of gas fractions among high Xsfr galaxies shows 
that feedback from starbursts does not necessarily result 
in a reduced Hi content. 

The Hi-richness of the starbursts indicates that they 
are not able to completely ionize or consume their neu¬ 
tral gas. Radiative feedback in high Xsfr galaxies does 
not appear to significantly affect their global HI gas frac¬ 
tions, perhaps because most of the HI gas mass resides 
at radii well outside the star-forming regions. The gas 
column near the starburst may be thick enough to pre¬ 
vent the escape of ionizing photons, and in addition, any 
decrease in HI from photoionization may be offset by an 
increase in Hi from H2 photodissociation. 


3.1.2. The Connection Between Hi Gas Fraction and 
Specific Star Formation Rate 

The ALFALFA Ha galaxies indicate that H1 gas frac¬ 
tion generally increases with sSFR, although the data 
show substantial scatter and this trend may flatten in 
the starburst regime. The relatively weak connection 
between sSFR or Ha EW and Hi gas fraction (Figure[8|) 
is seemingly at odds with the well-established tight cor¬ 
relation between Hi gas fraction and NU V-r color (e.g., 
ICatinella et al.l [20101 : Huang et alJ 120121b For instance, 
the GASS sample of high-mass galaxies shows that the 
best predictors of a galaxy’s Hi gas fra ction are NUV- 
r color and stellar mass surface density (ICatinella et alJ 
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Fig. 5. — A comparison of the ALFALFA Ha SFRs with SFRs corrected for extinction using (a) WISE 12 luminosities and ( b ) 22 /j, m 
luminosities. Gray circles indicate measurements using WISE point source photometry, and black crosses indicate galaxies with elliptical 
aperture photometry. The solid line shows a one-to-one relation. Representative error bars are shown in the lower right corner. 



Fig. 6.— Mhi/M* and M* for the ALFALFA Ha Fall sample. Symbols are described in Figure [4] The solid line shows a least-squares 
fit to the non-starburst galaxies. Two high-EW starbursts, shown in red, have Mhi/M* values 2 a higher than the best-fit line. These 
starbursts are discussed further in § 13.31 
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Fig. 7.— (a) Mhi/Lr vs. Esfr 5 with symbols as in Figure[4] Galaxies with higher Esfr have lower Mhi/Lr. ( b ) When we substitute 
M* estimates from SDSS photometry for Lr, Mhi/M* shows a weaker negative correlation with Esfr 5 and high Ssfr galaxies span a 
wide range of Mhi/M*. Color shows M*. 
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I20I0I) . Similarly, we observe a tight trend between HI gas 
fraction and SDSS g — r color (Figure 0i). The weaker 
trend with Ha-derived sSFR and strong trend with g — r 
suggest two possible explanations: ( 1 ) dust extinction 
drives the close correlation between g — r color and HI or 
( 2 ) galaxies’ global HI content is more closely linked with 
the SFR averaged over long ti mescales than with inst an- 
taneous star formation (e.g., iKannappan et all 1 201. '111 . 

To estimate the effect of dust extinction on SDSS color, 
we use the WISE 12 /im and 22 /jm luminosities (see 
to derive the E(g — r) color excess. Following 
IWen et all (IM1 . we calculate E(B — V) using 


vL v 10 0.4fc HQ xE(B-V) _ x 

^Ha(obs) ^ 


(4) 


where vL v is the 12 gm or 22 gm luminosity, L Ha (obs) is 
the Ha luminosity with no correction for internal ex¬ 
tinction, fee is the reddening curve at Ha, and a is 
a coefficient that depends on the extinction law and 
WISE band used (jWen et al.ll2014D . We t hen convert 
E(B — V ) to E(g — r), using the conversion in I Yuan et al.1 
(2013). The median Ay va lues inferred f o r the AL¬ 
FALFA Ha sample using the iCalzetti et al.1 (1200011 and 
iCardelli et, al.l (|1989f) extinction laws range from 0.45-0.6. 
These values are comparable to the face-on Ay values 
calc ulated from radiativ e transfer models of spiral galax¬ 
ies (iXilouris et all 119991: iBianchil l2007t iDe Gevter et al.l 
I2014H . In Figures 0> and 0, we show the relation be¬ 
tween Mhi/M* and extinction-corrected g — r color. Af¬ 
ter correcting the SDSS g — r colors for extinction, the 
trend with Hi gas fraction weakens substantially (Fig¬ 
ure 0 ■ 

Dust extinction has the most dramatic effect on the 
galaxy color of the re d, low sSFR galaxie s in o ur sam¬ 
ple. We therefore run iBruzual fe Chariot] (1200311 stellar 
population synthesis (SPS) models to confirm whether 
the above dust extinction corrections are realistic for 
low sSFR galaxies. We use the Pa dova 19 94 is ochrones 
at metallicities of 0.2 and 1 Zv,. dBressan et al.l 119931 : 
iFaeotto et al.lll 99~dll and the lChabrieil (12003H initial mass 
function. To estimate the expected colors of the low 
sSFR galaxies, we adopt a log-normal sta r form ation his¬ 
tory, as recommended bv lGladders et al.1 (2 013 for field 
galaxies. For present-day sSFRs of 10 -12 — 10” 11 yr -1 
and peak star formation at z = 0.7 — 3, the intrinsic 
g — r colors range from 0.5 to 0.8. While not exhaustive, 
these models hint at the range of colors we should expect 
for the reddest galaxies in the absence of dust. 

The SPS models cannot reproduce the reddest ob¬ 
served g — r colors, (g — r) > 0 . 8 , even assuming peak 
star formation at z > 6 . In fact, as suggested by the 
extinction-corrected colors in Figure 0 a wide range 
of star formation histories produces a relatively narrow 
range of stellar population colors. Colors redder than 
this range would then result from the effects of dust ex¬ 
tinction. The WISE extinction-corrected colors in Fig¬ 
ure El do appear more blue than predicted by the SPS 
models, with median g — r colors of 0 . 2 - 0 .3 using the 12 
/im band and 0.4-0.5 using the 22/mi band. Although 
WISE is less sensitive at 22/im th an at 12/im, the warm 
dust traced by the 22 /im band dJarrett et al.ll2013h may 
be a better indicator of the total stellar extinction than 
the polycyclic aromatic hydrocarbons (PAHs) traced by 


the 12/im band. Adopting a more realistic, bursty star 
formation history in the models could also lead to bet¬ 
ter agreement with the WISE extinction-corrected colors. 
Regardless, both the WISE data and the SPS models 
suggest that dust extinction could be responsible for the 
tight correlation between g — r color and HI gas fraction. 

Figures 0 and 0 show that, in star-forming galaxies, 
Hi gas fraction correlates weakly with sSFR; the tight 
connection between HI and g — r color may stem al¬ 
most entirely from a link between dust extinction and 
HI gas fraction. Nevertheless, we caution that ALFALFA 
misses most galaxies on the red sequence (jHuang et all 
i 2 im and including this population may result in a 
c learer trend betw een color and Hi gas fraction (cf. 
iGavazzi et al.ll20l3lf . 

For our sample, the correlation between Hi gas frac¬ 
tion and color may reflect an underlying trend between 
dust content and stellar mass. Figure [TU] demonstrates 
that the HYS Wderived ex tinc tions from the 12 /mi lu¬ 
minosities and ICardelli et al.1 (Il989fl extinction law cor¬ 
relate with stellar mass and anti-correlate with Hi gas 
fraction. Th e results f or the WISE 22 /im data and for 
the ICalzetti et all (12000( 1 extinction law are similar. The 
outlier at the highest E(B — V) in Figure [10] is UGC 
1488, which has a high inclination, prominent dust lane, 
and signs of morphological disturbance. Higher mass, 
star-forming galaxies have both higher dust masses and 
higher dust-to-gas ratios, due to the galaxies’ higher av¬ 
erage metallicities. Conversely, these same high-mass, 
dusty galaxies have low Hi gas fractions, due to higher 
H 2 conversion and star formation efficiencies or lower gas 
accret ion rates relative to their gas consumptio n rates 
(e.g., iDave et all 1 2 011 at IKannappan et al.ll2013H . 

The metallicities of galaxies set their dust content and 
ultimately drive the observed trend betwe en g — r color 
and H 1 gas fraction. Bothwell et al . (1201 3(1 show that the 
mass-metallicity relation of galaxies depends strongly on 
Hi content. As a galaxy evolves, the balance of Hi in¬ 
flows, star formation, and outflows determines both its 
metallicity and its Hi gas fraction. The tight trend be¬ 
tween color and H 1 gas fraction is therefore a manifesta¬ 
tion of the key role of H 1 gas flows in regulating galaxy 
metallicities. 

3.1.3. Kennicutt-Schmidt Law 

Although Hi gas fraction and sSFR do not appear 
strongly correlated, gas surface density may be the more 
relevant parameter for star formation . According to 
the Kennicutt-Schmidt Law (|Kennicuttl11998b ). galaxies 
with a higher SFR surface density should have a higher 
gas surface density. We calculate Esfr from Equation 0 
Due to the 3.5' resolution of the Arecibo beam, we do 
not have radii for the Hi gas. However, Hi diame¬ 
ters are ob served to scale with g alaxy optical diame ¬ 
ters (e.g.. IBroeils fc Rheel 119971: iSwaters et al.1 12002(1 . 
IBroeils fc Rheel (|1997H find that Hi radii are typically 
1.7 times larger than optical radii. We therefore use the 
SDSS radius containing 90% of the r-band light, Rqq to 
estimate the H 1 radius and calculate Sri as 

HI 27r(1.7i?9o) 2 " ( ) 

Since galaxies show substa ntial scatter in the ra tios of 
their Hi and optical radii (IBroeils fe Rheel fl997T) . these 
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Fig. 8. — (a)Mm/M* vs. sSFR, with symbols as in Figure |4] To compare with Figure [9b, we show ALFALFA Ha galaxies detected in 
WISE in black and non-detections in gray. (b) Mm/ M* vs. Ha EW. 


Ehi values are only rough estimates. 

We plot Ssfr and Shi for the sample in Figure fill We 
find no clear correlation between the Hi and SFR sur¬ 
face densities, as indicated by the low Spearman’s rank 
correlation coefficient of 0.16. This lack of correlation is 
consistent with results that show that molecul ar gas pri¬ 
maril y sets the Kennicutt-Schmidt Law (Schruba et al. 
l20TI . Most galaxies lie slightly below the threshold den¬ 
sity of ~10 M 0 pc~ 2 for conversion to molecular gas. On 
average, the starbursts tend to have elevated Shi relative 
to the rest of the sample; the high EW starbursts have 
an average Ehi more than 3 times higher than the non- 
starbursts, and the high sSFR starburst starbursts have 
an average Ehi 2 times higher than the non-starbursts. 
Higher surface densities may aid the conversion of atomic 
gas to molecular gas and fuel star formation. Alterna¬ 
tively, the starbursts may have more spatially extended 
H 1 than the other galaxies, or they may have lower 
metallicit ies a nd higher Hi saturation thresholds (e.g., 
IBolatto et al.l[20lTh . As with Hi gas fraction, the star- 
bursts’ Ehi values tend to be high, but they do not differ 
dramatically from the non-starbursts in the sample. 

Whether we consider sSFR, Ha EW, extinction- 
corrected g — r color, or Esfr, the Hi content of galax¬ 
ies shows only a weak connection with star formation. 
The starbursts have high H 1 gas fractions compared to 
the full ALFALFA Ha sample but show little to no in¬ 
crease in H 1 gas fraction relative to similar mass galaxies. 
These observations suggest that the Hi content of gas- 
rich galaxies remains relatively constant, even during a 
starburst episode. Any excess H 1 is efficiently converted 
into H 2 , and photodissociation of H 2 may balance the 
consumption or ionization of HI. 

3.2. Star Formation Efficiency 

Since the starbursts have slightly high, but not 
unusual, H 1 content for their masses, their current 
SFRs i mply a high H i-to-fF c onversion efficiency. 
Followin dSchiminovich et al.l (2 01(1 ), we refer to the SFR 
per unit H 1 gas mass as the H 1-based star formation ef¬ 
ficiency (SFE). This parameter indicates how efficiently 
galaxies are able to tap into their H 1 reservoirs and con¬ 
vert their Hi to molecular gas and stars. We calculate 


the inverse of the SFE, the Hi gas depletion timescale, 
as 


_ Mhi 
Ltep — SFR’ 


( 6 ) 


and we show the relation between t<jep and sSFR in Fig- 
ure[T2] The starburst H 1 depletion times are listed in Ta¬ 
ble □ Most of the sample shows no correlation between 
sSFR or EW and H 1 depletion time. However, the star- 
bursts tend to have short H 1 depletion times, despite the 
fact that they are Hi-rich. All the high EW starbursts 
and four out of six high sSFR starbursts have H 1 gas frac¬ 
tions above the sample median, and only two starbursts 
have depletion times below the sample median. In par¬ 
ticular, the two strongest or youngest starbursts, as mea¬ 
sured by EW, have both the highest H 1 gas fract ions and 
among the highest SFE. lSaintonge et al.l (1201 lbT) suggest 
that weak starbursts in the high-mass GASS sample may 
not be able to access H 1 easily, since they have the same 
H 1 fdep as non-starbursts. In contrast, our starbursts 
do show lower td ep , which may indicate that lower mass 
starbursts are efficiently converting H 1 to H 2 and stars. 

In general, td ep decreases slightly with M* (Fig¬ 
ure [l3k] i formasses below 10 10 M 0 (e.g., iBothwell et al.l 
l2009t IHuang et all 2012 1. The high H 1 gas fractions 
and low SFE observed for the lowest mass non-starburst 
galaxies may imply that their accretion rates exceed 
their rates of gas consu mption (e.g., IHuang et al.ll2012t 
Kannapp an et al.l I2013T) . either due to inefficient inward 
gas transport or inefficient H i-to-H 2 conversion. How¬ 
ever, the low-mass end also contains several galaxies 
with tdep an order of magnitude lower than the aver¬ 
age for their stellar mass. In addition to their higher 
scatter in tdep, the lowest mass galaxies also appear 
to ex hibit the largest scatter in sSFR (Figure (T3bb 
iLee et ahl (|2007tl explain the scatter in sSFR as an in¬ 
crease in“burstiness” at low stellar masses. 

We examine this change in scatter quantitatively in 
Figure [I3l In Figures fTfili and(T3b. we fit the mean tdep 
and mean sSFR as a function of M* using t he LOES S 
non-p aram etric local regressi on technique dClevelandl 
1197911 in R (|R Core Team [2014(1 . The LOESS method de¬ 
termines the mean at each M* by weighting data points 
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Fig. 9.— (a) SDSS g — r color and Mhi/M* show a tight correlation. Symbols and colors are the same as in Figure|8] The rr-axis is flipped 
to match the orientation of Figure|8] (6) Mhi/M* v s. SPSS g — r c olor, a fter correcting for extinction using the WISE 12 ^m luminosities. 
Bla ck points indic a te ext inctions derived assuming a lCalzetti et al.l (120001) extinction law, and blue points indicate extinctions derived with 
the ICardelli et al.l (119891) extinction law. The calculated errors include photometric uncertainties, the uncertainties in the coefficients in 
Equation [4| due to adopting different extinction laws, and the uncertainty in the E(B — V) to E(g — r) conversion, (c) Mhi/M* vs. SDSS 
g — r color, after correcting for extinction using the WISE 22 ^m luminosities. 


based on their distance from the M* value of interest. 
To illustrate the scatter about this mean, we plot the 
absolute value of the residuals from the mean fit in Fig¬ 
ures fT3fc andflOH. A LOESS fit to the residuals supports 
our rough by-eye assessment that the scatter about the 
mean in both tdep and sSFR is highest among the lowest 
mass galaxies. The gas-richness of star-forming dwarf 
galaxies, combined with their wide-ranging SFE, suggest 
that HI gas may accumulate until a dynamical distur¬ 
bance triggers gas inflows and compression. 

Galaxy structure may also play a role in determining 
Hi to H -2 conversio n efficiency and setting the Hi tdep- 
iBlitz fe Rosolowskvl (120061 ) argue that in disk galaxies, 
the midplane pressure sets the H 2 /HI ratio and hence the 
SFE. Since the midplane pressure shoul d sca le with stel¬ 
lar mass surface density, £* (e.g., IBlitz & Rosolowskvl 
12001 . we plot tdep as a function of £* in Figure [TT] We 


calculate 


where i?so,r is the SDSS r-band iPetrosianl (j!976T) half- 
light radius. Figure [14] shows that for the AL¬ 
FALFA Ha galaxies, td ep decreases with £*, as expected 
if higher surface density disks convert HI to stars more 
efficiently. 

Howev er, this trend is the exact opposite of that 
found bv lSaintonge et all (l2012f) for high-mass, optically- 
selected GASS galaxies, where the mean Hi td ep in¬ 
creases with £ ». We show the Hi td B n and of the 
GASS sample (jCatinella et al.l l2010. 20121 l2013l h calcu¬ 
lated following Equations [G] and 0 in Figure QjJ We 
only include GASS galaxies with Hi detections, and 
we use the “representati ve” GASS sample, w hich in¬ 
cludes gas-rich galaxies dCatinella et alJ 120131) . The 
GASS stellar masses are calculated from SDSS photom- 
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Fig. 10.— Mhi/M* vs. E(B — V) derived using the WISE 12 /mi luminosities and ICardelli et aD ^ 19891 ) extinction law. Color indicates 
stellar mass. 
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Fig. 11.— Ssfr vs. Shi 5 calculated from Equations [3] and [5] The values of Ehi are rough estimates, obtained by scaling the observed 
optical radii. The error bars account for statistical uncertainties only. 


etry (jCatinella et al.l I2013T) . We calculate SFRs from 
the tabulated GALEX NU V photometr y as descr i bed in 
ISchiminovich et aTI (I2010f) . assuming a iChabriei] (l2003h 
IMF. Among massive galaxies, UV -based SFRs agre e 
well with Ha-based estimates (e.g., iSalim et~ahl 120071) . 
which therefore permits a direct comparison with our 
sample. 

The discrepancy between the ALFALFA Ha galax¬ 
ies and the GASS galaxies is likely due to the quench¬ 
ing of star formation in more massive galaxies. Fig¬ 
ure M shows that the scatter in t<jep increases dra¬ 
matically above a £* of IQ 8 7 Mq kpc~ 2 , identifi ed as 
a “quenching threshold ” by iCatinella et al.l (1201011 and 
ISaintonge et al.l (|2011af) . The Hi depletion times of the 
few ALFALFA Ha galaxies in this high stellar surface 


density regi me are consistent wit h those of the GASS 
population. ISaintonge et al.l (12012H suggest that galaxies 
with high stellar surface densities may be more stable to 
fragmentat ion due to their bulge-do minated morpholo¬ 
gies (e.g., lOstriker fe Peebles!Il973h or may no longer 
have access to their HI reservoir. 

We use th e SDSS surface brightness profile fits of 
ISimard et al.l (|2011ll to quantify the morphologies of the 
GASS galaxies. These fits are not available for the 
ALFALFA Ha galaxies, which fall in the more recent , 
SDSS DR9 sky coverage. Following ISimard et al.l (1201 If) . 
we adopt a bulge+disk model if the probability that 
the bulge+disk model is not required is below 0.32; if 
the probability is higher, we use a pure Sersic model. 
Figure [15] shows that the GASS disk galaxies continue 
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Log Ha EW (A ) 



Fig. 12.— The relation between and Ha EW (left panel) and t<j ep and sSFR (right panel). The right panel only includes galaxies 
with reliable SDSS photometry. The dashed line in both panels indicates the Hubble time. Despite their high Hi content, most starbursts 
have shorter t^ep than the rest of the sample. 






Fig. 13.— (a) idep and (6) sSFR as a function of M*. The red line shows the mean value at each M*, as determined by a LOESS 
non-parametric local regression fit. The black lines show the 95% confidence intervals of the mean. In panels (c) and (d) we show the 
absolute value of the residuals about the mean fits in the upper panels. The red line shows a LOESS fit to the residuals, and the black 
lines again show the 95% confidence intervals. Low mass galaxies show higher scatter about the mean for both t^ ep and sSFR. 
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Log (M 0 kpc" 2 ) 


Fig. 14.— t ^ ep as a function of £* for the ALFALFA Ha sample (black circles, diamonds, and triangles) and the GASS sample (gray 
crosses). Below the “quenching threshold” of £* = 10 8 ' 7 M 0 kpc —2 <|Catinella et al.l 120101 : [Saintonge et alJl2011al) . t ^ ep decreases with 
increasing £*. The black cross indicates the median statistical errors for the ALFALFA Ha sample, and the dashed line indicates the 
Hubble time. 
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the trend of decreasing t<jep with S* seen in the AL¬ 
FALFA Ha sample, while spheroid-dominated galaxies 
deviate to higher tdep- Although all the offset galaxies 
appear to be bulge-dominated, they show no discernible 
trend between bulge-to-disk ratio and fdep- 

The spheroidal galaxies also appear offset to lower sS- 
FRs (Figure [16jj) compared to the rest of the GASS and 
ALFALFA Ha galaxies, supporting the idea that they 
have quenched their star formation. Interestingly, at a 
given £*, the Mhi/M* values of spheroidal galaxies are 
independent of sSFR (Figure H^bL indicating that the 
dispersion in fdep at the high £* end is not the result of 
a wide range in HI content. Instead, the low sSFR galax¬ 
ies appear unable to use their HI efficiently. In terms of 
absolute Mhi , galaxies with the highest HI masses tend 
to live just below the quenching threshold (Figure [Ht); 
galaxies with higher £* must have lost or consumed their 
HI gas during the quenching process. The lack of a trend 
with bulge-to-disk ratio or Mhi/-M* among the quenched 
galaxies suggests that differences in the spatial distribu¬ 
tion of the HI gas may account for their variation in 
SFE. 

While the spheroidal galaxies exhibit the highest 
fdep for their £*, starbursts have low fd ep relative 
to galaxies of a similar M* or £* (Figures HdlflGl) . 
Many studies have pointed out the role mergers may 
play in enhancing SFE (e.g.. iSolomon fe Sagel 119881: 
Combes et, alJll9f)4lBouche et alJl2007tlDi Matteo et all 

20071: iBournaud et all 1201 il l. Turbulence and gas flows 

during major mergers can lead to the efficient forma¬ 
tion of molecular ga s, shifting the ISM ga s distribution 
to hi gher densities (IBournaud et al.1 120111 : iPowell et al.l 
120131) . During a gas-rich merger, the total gas fraction 
of a galaxy may increase due to the addition of new gas; 
however, if the distribution of gas becomes skewed to 
higher gas densities, the H 2 gas fraction may increase 
more strongly than the H 1 gas fraction. This scenario 
could potentially explain the moderately high Hi gas 
fractions and short fdep of the ALFALFA Ha starbursts. 


3.3. Morphology and Mergers 

The starbursts’ morphologies and kinematics may 
demonstrate whether dynamical disturbances are en¬ 
hancing the SFE. One measure of morphological distur¬ 
bance is the 180° rotational asymmetr y. We calc ulate 
the /7.- band asymmetry, An, following iConselice et, al.1 
(I2000hlb as 


„ • l-fo - Jl8oK ,E \ B 0 ~ -B180K fo\ 

4 = mm(- -) — mm(- -), (8) 


EI J o 


EI J ol 


where Iq is the original image, /iso is the image rotated 
by 180°, B 0 is a background region, and /?i 8 o is the back¬ 
ground region rotated by 180°. We sum over all the im¬ 
age pixels within the Petrosian radius (lPetrosianll976i) at 
the i~i = 0.2 surface brightness level, where p is the ratio 
of the intensity at a given radius to the average intensity 
within that radius. We also scale the background region 
area to this radius. The minimization finds the rotational 
center that produces the lowest asymmetry value. Prior 
to the calculation, all images are background-subtracted, 
and we mask all stars within approximately five galaxy 


ra dii using the IRAlfl tas k imedit. As recommended 
by IConselice et al.1 (i2000bl h we adopt a signal-to-noise 
(S/N) cut of 100. Below this value, the scatter in An 
dramatically increases, and 41% of these low S/N AL¬ 
FALFA Ha galaxies have unphysical negative An values. 

We show images of the eight high EW and six high 
sSFR ALFALFA Ha starbursts and their calculated 
asymmetries in Figure [17] and list the asymmetry val¬ 
ues in Table [0 The eight starbursts with S/N > 100 
have asymmetries ranging from 0.18-0 .68, higher than 
the median sample asymmetry of 0.14. IConselice! (|2003f ) 
suggests that asymmetries above ~0.35 indicate major 
mergers. Figure [TH] shows the asymmetries of the sample 
as a function of Ha EW. All the starbursts have slightly 
elevated asymmetries, and two (ACC 330517 and AGC 
330500) are clearly major mergers. However, we find 
equally disturbed galaxies at lower values of Ha EW, 
indicating that not all H i-rich mergers are starbursts. 

If some of the starbursts are merging systems, they 
appear to be near coalescence. Simulations of major 
mergers also indicate that the peak s tar formation activ¬ 
ity should occur at th i s time (e.g . . IMihos fe Hernauistl 
11994 ICox et all 120081 : ILotz et ah I2010bl) . Intriguingly, 
the two most asymmetric starbursts, AGC 330517 and 
AGC 330500, are also the high EW starbursts with the 
highest offsets from the HI gas fraction and stellar mass 
relation (see Figure [6]). High gas-fraction mergers expe¬ 
rience more extended star formation, which c auses them 
to app ear asymmetric over a longer timescale (ILotz et al.1 
I2010af) . Alternatively, these two starbursts may be in an 
earlier merger stage than the other starbursts, prior to 
coalescence, when they have not yet consumed or ex¬ 
pelled much of their HI gas. The individual merging 
galaxies in AGC 330517 and AGC 330500 are still dis¬ 
tinguishable in Figure[T7l which supports this possibility. 
The other six starbursts with asymmetry measurements 
have An ~ 0.2 —0.3, only slightly higher than the average 
asymmetry of the full sample. These lower asymmetries 
could sim ply reflect clumpy star formation in these ga lax- 
ies (e.g., IConselice et aljfeoOObl : iFossati et akll2013fl . as 
appears to be the case for UGC 12821, or they could 
indicate a mild dynamical disturbance. The three most 
massive high sSFR starbursts, AGC 330186, UGC 470, 
and UGC 12821, show spiral structure and no signs of 
a major interaction. However, the presence of a nearby 
bright star blocks a full view of AGC 330186, and the 
galaxy does appear somewhat asymmetric in its struc¬ 
ture. Minor disturbances may cause the enhanced star 
formation in these spirals. AGC 120193, AGC 331191, 
and AGC 102643 show more obvious signs of morpholog¬ 
ical disturbance, with extended, asymmetric tails visible 
in their broadband images. The lower asymmetries in 
some of the starbursts could also b e consistent with a 
later merger stage; the simulations of lLotz et al.l (2 010 b) 
show that the peak star formation in major mergers may 
occur a few hundred Myr after the peak asymmetry. Mild 
asymmetries of An ~ 0.2 correspond to this peak star¬ 
forming, coalescence stage in mergers. If these less asym¬ 
metric high EW starbursts are indeed late-stage mergers, 

9 IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 





















































Hi and Star Formation in ALFALFA Ha 


17 





1.0 


0.8 


_o 

o 

0.6 * 

(n 

Q 

I 

O 

0.4 i 

D 

m 



0.2 


0 


Fig. 15.— t<j ep as a function of E» for the ALFALFA Ha non-starbursts (black circles) and starbursts (diamonds and triangles) and the 
GASS sample (colored circles and crosses). The GASS galaxies are colored by bulge-to-disk ratio from the lSimard et all (1 20111 ') SDSS fits. 
The crosses indicate GASS galaxies for which a pure Sersic profile is preferable; red crosses indicate a Sersic index > 2, and purple crosses 
indicate a Sersic index < 2. 


their lower Hi gas fractions relative to the more asym¬ 
metric starbursts AGC 330517 and AGC 330500 suggest 
that dynamical disturbances trigger the efficient conver¬ 
sion of HI into H 2 throughout the merger. The more 
asymmetric, earlier stage mergers, such as AGC 330517 
and AGC 330500, still exhibit excess Hi, which may be 
consumed or ionized by the end of the starburst phase. 

Disturbed gas kinematics may also be a sign of 
galaxy interactions. The Hi velocity profiles of merg¬ 
ing galaxies often exhibit asymmetries or wide, high- 
velocit y wings due to an excess of high-velocity ga s 
(e.g., (Gallagher et al.l fl981t iMirabel &; Sandersl I1988D . 
To quantify this high-velocity excess, we calculate the 
ratio of the HI velocity width at 20% of the peak flux to 
the width at 50% o f the p eak, W 20 /W 50 , as suggested 
by iConselice et al.1 (2000a). Since the ALFALFA ve¬ 
locity resolution is 10 kms -1 , we calculate an upper 
limit to W 20 /W 50 for galaxies with W 20/2 — W^q/2 <10 
kms -1 . We show the Hi width ratios and Ha EWs of 
the sample in Figure [T9] and Table [Tj All of the high EW 
starbursts have H 1 width ratios greater than the sam¬ 
ple median, although as with morphological asymmetry, 
high Hi width ratios occur in no n-starbursts as we ll. In 
simulations of merging galaxies. iPowell et al.l (120131) find 
that interaction-driven turbulence, rather than super¬ 
nova feedback, creates the high gas velocity dispersions 
observed in mergers. They argue that this enhanced tur¬ 
bulence increases the SFR by shifting the gas distribution 
to higher densities. The observed Hi gas disturbances 


are therefore a possible cause of the high star formation 
efficiencies in the high EW starbursts. 

Non-interacting dwarf galaxies may also show broad 
wings in their Hi profiles, however, and the high 
W 20 /W 50 ratios of the starb ursts do not necessarily prove 
they are interacting (e.g., Gallagh er et al.lH98lD . We 
display the Hi velocity profiles of the high EW and 
high sSFR starbursts in Figure l20l While the Hi pro¬ 
files of several starbursts appear broad or irregular, the 
profiles of some of the lower-mass starbursts, especially 
AGC 112546 and AGC 122866, appear quite narrow and 
only have upper limits on W 20 /W 50 . For comparison, we 
show examples of non-starburst Hi profiles in Figure Ell 
For each starburst, we randomly selected a non-starburst 
counterpart with M * within ^0.5 dex and with a similar 
S/N Hi spectrum. Since one starburst, AGC 120193, 
does not have SDSS photometry, we selected a non- 
starburst with Mhi within 0.5 dex. The Hi velocity pro¬ 
files of the four most massive high EW starbursts, AGC 
330517, AGC 120193, AGC 330500, and AGC 331191, 
do appear to have broader wings than the straight-sided 
profiles of the non-starbursts. However, the profiles of 
the four lower mass high EW starbursts (M* < 1O 8 M 0 ) 
are indistinguishable from non-starburst profiles of simi¬ 
lar mass galaxies. We caution that the resolution of the 
H 1 spectra may be too low to accurately probe the pro¬ 
file shape s of these galaxies. Nevertheless, as discussed 
in § 13.1.11 the low-mass starbursts also have H 1 gas frac¬ 
tions comparable to those of low-mass non-starbursts. 
Therefore, the global Hi gas fraction may not be the 
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Fig. 16. — tdep as a function of E* for the ALFALFA Ho; non-starbursts (circles) and starbursts (diamonds and triangles) and the GASS 
sample (crosses). Color indicates sSFR (top panel), Mhi/M* (center panel), and Mhi (bottom panel). 
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Fig. 17.— Morphologies of the eight high EW ALFALFA Ha starbursts (top panel), organized in order of Ha EW (see Table [TJ. The 
highest EW starburst is at the upper left, and the lowest EW starburst is at the lower right. The bottom panel shows the morphologies 
of the six high sSFR starbursts organized by decreasing sSFR. The ALFALFA Ha R-band images appear in red, SDSS g-band images 
appear in green, and the continuum-subtracted Ho; emission is shown in blue. All images have a logarithmic brightness scale. AGC 120193 
does not have SDSS images. The solid white line in each panel corresponds to 5 kpc. The R-band asymmetries are indicated for the eight 
starbursts with S/N> 100. The yellow region above AGC 330500 is a bright star. 
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Fig. 18.— The ij-band asymmetries and Ha EWs for galaxies with S/N> 100. The solid line shows the median asymmetry and the 
dotted line shows the Ar > 0.35 criterion for identifying major mergers. The downward arrow indicates an outlier with unphysical negative 
asymmetry. For clarity, two overlapping data points for the high sSFR starbursts have been artificially displaced in EW. 



Log Ha EW (A) 


Fig. 19.— Hi velocity width ratios, W 20 /IF 50 , and Ha EWs of the ALFALFA Ha sample. The dashed line shows the median W 20 /W 50 ', 
this median includes galaxies with upper limits on 11 20 /Hr,a, and therefore may be overestimated. The starbursts, shown by red diamonds 
and triangles, tend to have above average HI width ratios. Representative error bars are shown at the upper left of the plot. 


primary physical parameter that sets the level of star 
formation in low-mass galaxies. Infalling gas clouds are 
one p ossible t rigger o f star f ormation in dwarf gal axies 
(e.g; iGordon fc Gottesmanll 1 98 It iLo ncz-Sa nchez et al.1 
l2012t IVerbeke et al.l 1201411 . IVerbeke et al.l (12014 )' show 
that this type of interaction may generate multiple star- 
burst episodes; disturbed gas kinematics may only be 
observable during the initial starburst and with favor¬ 
able viewing orientations. These low-mass starbursts 
may therefore no longer show H 1 disturbances and may 
have consumed some of their Hi during a prior burst. 

The high sSFR starbursts generally show more ordered 
Hi velocity profiles than the high EW starbursts. As 
with the some of the high EW starbursts, two galax¬ 
ies (AGC 122420 and AGC 102643) have inconclusive 
Hi profiles due to low resolution or low S/N. The other 
high sSFR starbursts show no obvious disturbances. 
The more massive spirals AGC 330186 and UGC 470 


have straight-sided profiles indicative of orderly rotation, 
while the Gaussian profile o f UGC 12821 is_ likely due to 
its face-on orientation (e.g., lHavnes et al.lll998lf . Finally, 
the H i-rich, low surface brightness galaxy AGC 320466 
also shows ordered rotation in its Hi gas. With a high 
Mhi/M* = 25, much of AGC 320466’s Hi gas may not 
reside near its star-forming regions. Unlike the broad 
profiles of several of the high EW starbursts, the ordered 
velocity profiles of most of the high sSFR starbursts do 
not indicate any major disturbances. 

Although the traditional picture of merger-driven star 
formation emphasizes nuclear starbursts, recent studies 
have suggested that mergers may also cause a su b stantial 
increase in extended sta r formation (e.g., lEllison et al.1 
l2013l : IPowell et ahll2013l) . In Figure [251 we show the frac¬ 
tion of star formation within R 50 for starbursts and non- 
starbursts. We only include galaxies with a S/N> 10 
in the Ha photometry. Three starbursts (AGC 112546, 
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Fig. 20.— The Hi velocity profiles of the high EW starbursts (a) and high sSFR starbursts (h). ordered as in Figure 1171 


AGC 122866, and AGC 331191) show highly concen¬ 
trated star formation, with more than 80% of their 
Ha emission contained within f? 5 o- However, other star- 
bursts, particularly AGC 120193 and AGC 333529, show 
H II regions offset from the main R-band center of the 
galaxy. The majority of the star formation in these two 
starbursts is outside R 50 . Gas flows in starbursts fuel 
nuclear star formation in many, but not all, cases. 

4. DISCUSSION 
4.1. Hi and Star Formation 

The ALFALFA Ha galaxies confirm previous obser¬ 
vations that at a given stellar mass, more H 1-rich 
galaxies typically have higher SFRs ()Wang et al.l 1201 it 
iHuang et alJl2012 ). Since a high Hi gas fraction relative 
to other galaxies of a similar mass may in dicate th e pres - 
ence of recently accreted gas (e.g., iMoran et al.l l2012h . 
this trend may suggest that Hi accretion actively fuels 
star formation in low-redshift galaxies. Nevertheless, the 
scatter in the Mhi/SFR ratio at a constant stellar mass 


shows that additional factors affect the ability of galaxies 
to access fresh H 1 fuel and efficiently convert it to H 2 . In 
particular, the starbursts generally show high H 1-based 
star formation efficiencies (i.e., short H 1 depletion times) 
compared to similar mass galaxies (§ liOll . We also find 
that the starbursts’ Hi gas fractions are not unusually 
elevated for their stellar mass, despite their much higher 
SFRs IS IS. 1. ID . Therefore, in most cases, a high efficiency 
of H 1 to H 2 conversion, rather than a large influx of gas, 
is generating their high levels of star formation. Dy¬ 
namical disturbances are a likely cause of this enhanced 
efficiency, allowing the starbursts to transport Hi gas 
inward and generate high gas column densit ies m or e eas¬ 
ily than non-starbursts (e.g., iBarnes fe Hernauistlll99ll : 
iDi Matteo et al.lfeOOTt iHopkins et al.ll2013l) . 

The balance between H 1 infall and feedback may lead 
to a relatively stable H 1 supply in gas-rich starbursts. 
Addition of H 1 gas through interaction with a gas cloud 
or gas-rich galaxy may initially raise the H 1 gas fraction. 
This excess H 1 quickly disappears, however, as gas com- 
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Fig. 21.— The H I ve locity profiles of 14 non-starbursts, selected to have similar stellar masses and S/N spectra as the corresponding 
starbursts in Figure l20l Since one starburst (AGC 120193) does not have a stellar mass estimate, we select the corresponding galaxy (UGC 
1478) on the basis of HI mass. 


pression or turbulence resulting from these interactions 
efficiently converts the Hi gas to H 2 and stars. As the 
starburst progresses, the Hi gas fraction may stabilize 
instead of continuing its decrease. The starburst may 
not affect Hi gas in the outer galaxy, and photodissoci¬ 
ation of molecular gas may replenish H 1 in the inner re¬ 
gions. Indeed, previous observations have demonstrated 
that photodissociated gas may be a significant compo¬ 
nent o f the ISM in the inn er regions of starburst galaxies 
(e.g., iStacev et all 119911) . Radiative feedback does not 
appear sufficient to deplete the H 1 gas in the starbursts 
through ionization, and the two high EW starbursts with 
the highest ratios of Ha luminosity to H 1 mass (i.e., the 
shortest tdep) also show the highest Hi gas fractions. Af¬ 
ter the starburst, the remaining Hi reservoir may con¬ 
tinue to fuel star formation, or additional feedback mech¬ 
anisms, such as supernovae, may ultimately quench the 
burst. 


4.2. The Link Between Hi, Metallicity, and Dust 

For the ALFALFA Ha sample as a whole, we find 
that galaxies’ sSFRs and Hi gas fractions are only 
weakly correlated (§ 13.1.111 . This result contrasts with 
the strong observed correlation between galaxy color 
and Hi gas fraction identified in previ o us studies (e.g., 
Kannappanl |2004 IZhang et al.l l2009t iCatinella et alJ 
2010HHuang et al1l2012D which suggests th at galaxy sS¬ 
FRs are linked to their Hi content (e.g., IZhang et all 
120091 : iHuang et, al.l l2012t) . We likewise observe a tight 
correlation between H 1 gas fraction and g — r color. For 
our sample of gas-rich galaxies, we demonstrate that dust 
extinction, rather than star formation, drives the tight 
trend between H 1 and color. 

In essence, the Hi and color relation is a manifesta¬ 
tion o f the well-known galaxy mass-metallicity re lation 
(e.g., ILeaueux et al.lH979l : iTremonti et all 120041) . Re¬ 
cent studies suggest that galaxies in fact lie on a fun¬ 
damental plane of stellar mass, SFR, and metallicity 
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Fig. 22.— The fraction of Ha emission within the R-band half-light radius for non-starbursts (solid lines) and starbursts (dashed lines). 
Only galaxies with an Ha S/N> 10 are shown. 


(iLara-Lonez et al.ll2010tlManmicci et al.ll20l6h . known as 
the ‘fundamental metallicit y relation’. Gas content likely 
drive s this relation (e.g., jPave et al.l 120121 iLillv et all 
I2013H . and iBothwell et al.1 ( 2013T) show that a funda¬ 
mental relation exists between Hi mass, stellar mass, 
and metallicity. In this model, Mhi/M* and metallic¬ 
ity should anti-correlate. Since more metal-rich galax¬ 
ies have higher dust content, we would therefore expect 
galaxies with high Mhi/M* to have less dust extinction 
and bluer colors. 

Several processes could produce the proposed relation 
between galaxy stellar mass, Hi content, and metallicity. 
Over their lifetimes, galaxies with higher stellar masses 
will have both consumed more HI gas and produced more 
metals via star formation. In addition, the higher poten¬ 
tial wells of more massive galaxies should allow them to 
retain metals more effectively (e.g., iDave et akll2011al) . 
At a given stellar mass, galaxies that have experienced a 
recent inflow of metal-poor gas from the IGM should have 
lower overall metallicities as well as higher Hi content. 
Finally, we note that the elevated dust content of metal- 
rich galaxies will also allow them to convert their HI gas 
to H 2 more efficien tly by shielding molecular gas from 
UV radiation (e.g., IKrumholz et all (20091 iBolatto et, al.1 
HolO), thereby reducing their Hi fractions further. This 
connection between H 1 content and dust extinction pro¬ 
vides a natural explanation for the observed relation be¬ 
tween H 1 gas fraction and color. H 1 inflows and outflows 
are fundamentally important to the evolution of galaxy 
metallicities, while galaxy metallicity may also influence 
the H 1 supply by promoting H 2 formation. 


4.3. Galaxy Structure and Hi Conversion 

Another factor that may enhance the conversion of H 1 
to H 2 is galaxy structure. We find that the Hi deple¬ 
tion time in disk-dominated systems anti-correlates with 
stellar surface density (§ HOI) . This relation may indi¬ 
cate that higher disk mid-plane pressure aids the for¬ 
mation of molecular clo uds a nd in creases the H 2 /H 1 ra¬ 
tio, as proposed by iBlitz fe Rosolowskvl (l2006lh How¬ 
ever, alternate t heories of H? formati on based on self- 
shielding (e.g., IKrumholz et al.l l2009l i may re sult in a 
simil ar stellar surface density scaling relation (|Fu et al.1 
[MS). Early supernova feedback in low-mass galaxies 
may lead to the delayed infall of Hi gas, resulting in an 
increased H 1 content and l onger depletion time at the 
present day dFu et all 12010 1. By suppressing past star 
formation episodes, this efficient feedback should also 
lead to lower present-day stellar surface densities. 

We only observe this correlation between HI depletion 
time and stellar surface density in disk galaxies, how¬ 
ever. We show that high- mass, early-ty pe galaxies from 
the GASS sample (ICatinella et al.ll2013D do not follow a 
similar relation. Within this early-type sample, we find 
no correlation between H 1 depletion time and morphol¬ 
ogy, which suggests that bulge strength is not affecting 
the current SFE in these galaxies. Instead, we suggest 
that the scatter in SFE results from differences in the 
spatial distribution and surface density of the H 1 gas. 

4.4. Mergers and Hi Conversion 

For several of the starbursts, major mergers may in¬ 
duce their high SFRs and high HI to H 2 conversion ef- 
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ficiencies (TO- During interactions, tidal torques pro¬ 
pel gas toward the center of gal axies, t h ereby raisin g 
the gas column density (e.g., I Barnes fe Hernauistlll991lf 
and allowing galaxies to access their external Hi sup¬ 
ply. Furthermore, turbule nt motions during the merger 
may compress gas clouds (lElmegreen et al-l 11993 ). lead¬ 
ing to a higher fraction of dense gas and a higher 
H 2 /Hi ratio dPowell et al-l 120131) . Starbursts also ap¬ 
pear to form sta rs fr om H 2 more efficiently than non- 
starbursts (e.g., iKennicuttlIl998bh . perhaps due to the 
highe r mean density and sh orter free-fall time in their 
ISM (jKrumholz et all 12012! ). A merger scenario could 
therefore explain the moderate Hi but high SFRs of at 
least some of the ALFALFA Ha starbursts. 

Evidence of morphological and kinematical distur¬ 
bances support a merger origin for several of the star- 
bursts. We find that the high EW starbursts’ opti¬ 
cal morphologies are more asymmetric and their HI ve¬ 
locity profiles have wider wings than most of the 
non-starbursts. Merger simulations predict that the 
peak star formation act ivity should occur near coales¬ 
cence (e.g. , IMihos fe Hemciuistl 11994 ICox et all 120081 
ILotz et all I2010bh . consistent with the starbursts’ mor¬ 
phologies. The highest asymmetries, however, should ap¬ 
pear shortly before coalescence. Interestingly, our two 
most asymmetric starbursts, which appear to be in this 
merger stage, also have higher-than-average HI gas frac¬ 
tions for their stellar mass. Their high Hi gas frac¬ 
tions suggest that HI conversion to H 2 may be an on- 
goin g process dur ing me rger-d riven star formation (e.g., 

I Hibbard fc van Gorkoml[l996l) . 

Due to the unresolved nature of the Arecibo Hi ob¬ 
servations, we do not know the surface density of the 
Hi gas in the ALFALFA Ha starbursts. However, we 
find that they have above-average ratios of Mhi to the 
galaxy optical area. Therefore, either they lie closer to 
the maximum column density threshold for H I, which 
would aid H 2 formation, or their HI is significantly more 
extended than the optical disk. For instance, extended 
HI tidal tails may contain_a la rge fractio n of the HI gas 
in merging galaxies (e.g., iHibbard fc van Gorkomlll996D . 

The evidence for mergers among the lowest mass star- 
bursts is less clear, in part due to their lower S/N opti¬ 
cal images and Hi spectra. Their Hi velocity profiles 
do not noticeably differ from the profiles of low-mass 
non-starbursts, which may indicate a lack of kinemati¬ 
cal disturbances. Resolved HI observations likewise indi¬ 
cate that ord ered kinematics are not uncommon in dwarf 
starbursts dLelli et al-ll2014T) - However, if the dwarf star- 
bursts experience periodic bursts, jury initial disturb ance 
may no longer be evident (e.g., lYerbeke et al.ll20f4l ). We 
observe an increased scatter in SFE at the low-mass end 
of our sample, which supports this scenario of episodic 
bursts in dwarf galaxies. Supernova feedback may have a 
stronger effect in the low potential wells of dwarf galax¬ 
ies, leading to temporary quenching and a renewed burst 
of sta r formation as neutral gas falls back i nto the galaxy 
(e.g., iLee et alll2f)f)7t lYerbeke et al.ll2014fl . Highc r reso¬ 
lution imaging and HI spectral observations will be nec¬ 
essary to determine whether or not the dwarf starbursts 
are merging systems. 

Finally, the spiral structure and orderly HI veloc¬ 
ity profiles of the three highest mass starbursts, AGC 
330186, UGC 470, and UGC 12821, suggest they are not 


experiencing a major merger. UGC 470 has a higher than 
average gas fraction for its stellar mass and is known 
to have an unusually extended Hi disk (iDowelllfeOlOl ). 
It is also the only starburst in the sample with a typi¬ 
cal Hi depletion time for its stellar mass or Esfr (Fig¬ 
ures [13] and [J4|. Unlike the other starbursts, UGC 470’s 
excess gas supply, rather than a high Hi-to-H 2 conver¬ 
sion efficiency, may account for its enhanced star forma¬ 
tion. Minor mergers or other minor disturbances may 
temporarily increase the SFR for AGC 330186 and UGC 
12821, and most of the HI disk may remain undisturbed. 

4.5. Hi Cycles in Starbursts 

We therefore propose the following picture for the cy¬ 
cling of HI and star formation throughout the strongest 
starburst episodes. In the intermediate mass starbursts, 
a gas-rich major merger triggers the star formation. As 
the two galaxies approach coalescence, the total stellar 
mass of the system increases. Since the two individ¬ 
ual galaxies have HI gas fractions typical of lower mass 
galaxies, the combined system appears to have a higher 
Hi mass than other galaxies of a similar total mass. 
Strong tidal torques drive HI gas outward to form tidal 
tails and inward to fuel star formation. At this time, the 
system appears progressively more morphologically dis¬ 
turbed, and the increased gas flows and turbulence create 
kinematically disturbed HI profiles. The two most asym¬ 
metric starbursts in our sample, AGC 330517 and AGC 
330500, are possible examples of this stage, exhibiting 
disturbed morphologies, disturbed Hi kinematics, and 
elevated Hi gas fractions. As the Hi gas flows inward, 
turbulence and gas compression efficiently convert HI to 
H 2 , reducing the HI gas fraction and increasing the SFR. 
Although the gas flows bring Hi inward, star forma¬ 
tion is not necessarily restricted to the nuc lear reg ion, 
consi stent with recent simulations (e.g., IPowell et all 
120131) and with the varied Ha morphologies of the AL¬ 
FALFA Ha starbursts 1? 13.31) . Turbulent, dense clumps 
may arise throughout the merging disks, and the precise 
merger configuration affects the spatial distribution of 
star formation. The influx of HI gas raises the ISM col¬ 
umn density, and the star formation rate increases to the 
point where supernova-drive n turbulence supports t he 
enhanced weight of the ISM (lOstriker fe Shettvll201 lh . 

Peak star formation occurs near final coalescence, 
as the morphological disturbances are fading (e.g., 
ILotz et alJl2010br i. which may explain the slightly lower 
asymmetries of some of the starbursts (e.g., AGC 120193 
and AGC 331191). However, at the time of the peak 
SFR, the turbul ent moti ons d riving star formation are 
still high (IPowell et all I2013T) and the HI kinematics 
should still show higher velocities, as we observe in the 
high velocity wings of the high EW ALFALFA Ha star- 
bursts’ Hi profiles. The Hi content also begins to drop, 
due to the enhanced H 2 fraction and photoionization. 
Nevertheless, the starbursts still maintain a large HI sup¬ 
ply, comparable to similar mass non-starburst galaxies 
(§ 13.1.11) . Radiative feedback does not completely ion¬ 
ize the starbursts’ Hi reservoirs and may only ionize 
the H I gas near the starburst region (e.g., lHanish et all 
120101 ). In addition, H 2 photodissociation may compen¬ 
sate for the ionization of HI in the inner regions. Consis¬ 
tent with this scenario, recent HI observations of post¬ 
merger galaxies demonstrate that star formation and 
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feedback do not n oticea bly dep lete th e Hi reservoirs of 
merging galaxies (|Ellison et al.l 120151) . Ultimately, gas 
consumption or feedback terminates the starburst, and 
the final merged galaxy exhibits a higher stellar mass, 
higher metallicity, and lower HI gas fraction than its in¬ 
dividual progenitors. 

Lower mass galaxies may experience recurrent star- 
bursts, leading to the enhanced scatter in td ep and sSFR 
in the low-mass end of our sample. They may accu¬ 
mulate HI from the IGM until a dynamical disturbance 
initiates strong star formation. For instance, this initial 
trigger could be a merger with anothe r galaxy, an inter¬ 
action with a dark matter subhalo (jHelmi et al.l 120121) . 
or a merger with a gas cloud. In the case of a gas cloud 
merger, the initial interaction should raise the HI content 
and disturb the HI ki nematics, b ut may not trigger an 
immediate starburst (lYerbeke et alJf20l4 h Subsequent 
infall from the gas cloud or enhanced turbulence may 
then induce a later starburst. The morphology of the 
star formation should vary depending on the nature of 
the interaction and the gas cloud trajectory, con siste nt 
with the range of Ha morphologies we observe (§ & 

Regardless of the cause of the initial burst, supernova 
feedback rapidly quenches the initial starb u rst d ue to the 
low potential well of the galaxy (e.g.. iDekel fe SilldlTfJSfil 
IStinson et al.l 1200 7t IHopkins et al. 1201411 . This feedback 
may increase the HI velocity dispersion, but these kine- 
matic disturbances s hould lag the peak star formation 
(jVerbeke et al.l I20l4lh After the starburst dies down, 
the re-accretion of previously expelled gas may t rigger 
a series of subsequ ent starbursts (e.g., iLee et all 120071 : 
IStinson et al.l 120(371) . In addition, massive star clusters 
formed in a prior burst may gen erate tor q ues t hat con- 
tinue to drive HI gas inward (e.g., lElmegreen et al.ll2012T ) 
to fuel future starbursts. Thus, low-mass starbursts 
may lack the clear signs of kinematic disturbances that 
should characterize higher-mass interacting starbursts, 
as illustrated by the similar Hi velocity profiles of the 
low-mass starbursts and non-starbursts in our sample 
(§ & These starburst cycles may gradually reduce 
a dwarf galaxy’s HI reservoir but are unlikely to de¬ 
plete it entirely, given the starbursts’ high Hi gas frac¬ 
tions. As with the more massive starbursts, photoion¬ 
ization may not affect the outer HI gas and may be 
balanced by photodissociation near the starburst. Con¬ 
sequently, the Hi gas fractions of low-mass starbursts 
may also remain relatively constant, with little variation 
with respect to non-starbursts, and may provide a plen¬ 
tiful neutral gas supply capable of fueling multiple star- 
burst episodes. The similar Hi gas fractions of the AL¬ 
FALFA Ha starbursts and non-starbursts agree with this 
scenario IS 1.3.1.111 . as do observations of the Hi content 
of interacting dwarf galaxies relative to isolated dwarfs 
(iStierwalt et al.ll2014f ). By quickly quenching star forma¬ 
tion, the high feedback efficiencies in the dwarf starbursts 
may prevent rapid increases in metallicity or stellar mass. 
This scenario is consistent with the flat star formation 
histories and relatively i neffici ent sta r formation inferred 
for dwarf galaxies (e.g., iBehroozi et al.ll2013( ). 

The dwarf starbursts’ periodic star formation mode 
may be particularly relevant to high-redshift star forma¬ 
tion. Milky Way progenitors at z ~ 2 may accrete gas 
from the IGM i n discrete episodes that trigger enhanced 
star formation (IWoods et al I [20~i~il) . However, as with 


the dwarf starbursts, feedback efficientl y expe ls t he _gas 
and delays its consumption (e.g., IWoods et al.ll2014lf . 
As a result, these galaxies may have variable SFRs and 
may maintain an el evated Hi content t hat fuels later 
star formation (e.g., IHopkins et al.ll20T4 j. However, un¬ 
like z = 0 galaxies, z = 2 galaxies exper ie nce higher 
accretion rates from the IGM (e.g., iKeres et~akl 120051 ). 
High-redshift galaxies may therefore have higher average 
SFRs and larger gas reservoirs than we observe for the 
ALFALFA Ha sample. Feedback and variable SFRs in 
low-mass galaxies also have important impli c ations for 
the reionization of the IGM. IWvithe fc Loebl (1201,1( 1 ar¬ 
gue that by suppressing star formation, efficient feedback 
may reduce the contribution of the lowest mass galaxies 
to reionization. Constraining the starburst duty cycle in 
dwarf galaxies is therefore important to understand both 
galaxy star formation histories and reionization. 

5. SUMMARY 

The ALFALFA Ha survey represents the first opportu¬ 
nity to compare the properties of gas-rich starbursts and 
non-starbursts within a statistically uniform H i-selected 
sample. In this work, we analyze the HI gas fractions, 
Hi depletion times, Hi kinematics, and optical mor¬ 
phologies of 14 starbursts within the 565 galaxies that 
make up the ALFALFA Ha Fall-sky sample. This sam¬ 
ple illuminates the roles of gas accretion and feedback in 
determining the HI content of starburst galaxies and in 
triggering and sustaining star formation. 

Our main results are as follows: 

1. On average, the ALFALFA Ha galaxies with higher 
instantaneous sSFRs tend to have slightly higher 
HI gas fractions, but this trend is weak and shows 
substantial scatter. Galaxies with sSFRs that differ 
by an order of magnitude may still have the same 
HI gas fraction, and most of the starburst galaxies 
have HI gas fractions similar to galaxies with signif¬ 
icantly lower SFRs. In contrast, we observe a tight 
trend between HI gas fraction and g — r color. We 
show that dust extinction, rather than recent star 
formation, is primarily responsible for the tight HI- 
color correlation. This link between dust extinction 
and Hi gas fraction likely stems from the relation 
between stellar mass, metallicity, and HI content 
in galaxies (iBothwell et al.l[20111 ). 

2. Disk galaxies lie on a sequence of decreasing HI de¬ 
pletion time with increasing stellar surface den¬ 
sity. The observed trend is consistent with the 
idea that higher midplane pressures encourage 
the formation of H 2 from H 1 dBlitz fe Rosolowskvl 
120061). Disk galaxies from the GASS sample 
dCatinella et alJT2010t I2013H also fall on this se¬ 
quence, while spheroid-dominated systems are off¬ 
set to higher fdep ■ The spread in the HI depletion 
times of the spheroids reflects a spread in sSFR, 
but shows no trend with either Hi gas fraction or 
bulge-to-disk ratio. Instead, the spatial distribu¬ 
tion of H 1 in early-type galaxies likely determines 
whether gas clouds can reach the necessary densi¬ 
ties to form molecular gas. 

3. Gas-rich starbursts are able to maintain a relatively 
constant H 1 supply. Most of the 14 starbursts show 
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little to no increase in Mhi/M* or Ehi relative to 
galaxies of a similar mass; by the time of the star- 
burst episode, any excess atomic gas has already 
been converted into H 2 . However, we do find a few 
exceptions to this scenario. The extended HI disk 
of UGC 470 may fuel its elevated star formation. 
In addition, the two most optically asymmetric 
starbursts, which appear to be in a pre-coalescent 
merger stage, do show some evidence for enhanced 
Hi gas fractions. These asymmetric galaxies sug¬ 
gest that the conversion of excess Hi to H 2 is an 
ongoing process during mergers. 

4. Ionization does not appear to substantially de¬ 
plete the starbursts’ Hi gas, and photodissocia¬ 
tion of H 2 may compensate for decreases in H 1 due 
to consumption and ionization. Although we do 
not find that starbursts are unusually H i-rich, we 
also do not find that starbursts are H i-deficient, 
as suggeste d for starbursts in the SINGG sample 
(jOev et al.ll2007fl . Instead, high Esfr galaxies span 
a wide range of HI gas fractions. The similar¬ 
ity of the Hi gas fractions of starbursts and non- 
starbursts may indicate that the intense ionizing 
radiation of the starbursts does not penetrate to 
the outermost regions hosting much of the H 1 mass. 

5. The starbursts use their Hi more efficiently than 
the rest of the sample, as indicated by their lower 
H 1 depletion times relative to galaxies of a similar 
mass or stellar surface density. Major mergers may 
cause these high efficiencies in at least some of these 
starbursts, as suggested by the starbursts’ asym¬ 
metric optical morphologies and the wide H 1 veloc¬ 
ity profile wings in several starbursts. The high op¬ 
tical asymmetries of two starbursts are consistent 
with major mergers. The lower, but slightly ele¬ 
vated asymmetries of an additional six starbursts 
may indicate either more minor disturbances or 
mergers near coalescence, the merger phase pre¬ 
dicted to cause the largest enhancement in SFR. 
Finall y, consistent with recent merger simulations 
(e.g., iPowell et al.l l2013h . we find that extended, 
rather than nuclear, star formation may dominate 
the morphologies of some starbursts. 

6. While some of the starbursts are likely mergers, 
the lowest mass starbursts, with M* < 1O 8 M 0 , do 
not show clear evidence of disturbed optical mor- 
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phologies or HI kinematics. These dwarf starbursts 
may undergo periodic bursts, possibly triggered by 
a previous interaction, in which case un usua l gas 
kinem atics might not be apparent dVerbeke et al.1 
l2014j l. Large fluctuations in SF R appear charac ¬ 
teristic of dwarf galaxies (e.g., iLee et ali l2(Y)l . 
and the dwarf galaxies in our sample also show a 
larger scatter in SFE than more massive galaxies. 
The large, apparently sustainable H 1 gas fractions 
of low-mass starbursts may provide ample fuel for 
multiple generations of starbursts. 

The ALFALFA Ha galaxies demonstrate that while 
starbursts may differ dramatically from non-starbursts in 
their molecular gas content, the atomic gas fractions of 
starbursts and non-starbursts are similar. Efficient con¬ 
version of atomic to molecular gas reduces any H 1 excess, 
and the localized starburst may not affect the galaxy’s 
extended Hi reservoir. The Hi gas fractions of most 
low-mass starbursts may remain approximately constant 
throughout the burst, and only the most extreme star- 
bursts at low redshift may significantly disrupt their 
H 1 gas supply. 
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